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Abstract 

We describe how the entanglement renormalisation approach to topo- 
logical lattice systems leads to a general procedure for treating the whole 
spectrum of these models, in which the Hamiltonian is gradually simpli- 
fied along a parallel simplification of the connectivity of the lattice. We 
consider the case of Kitaev's quantum double models, both Abelian and 
non-Abelian, and we obtain a rederivation of the known map of the toric 
code to two Ising chains; we pay particular attention to the non-Abelian 
models and discuss their space of states on the torus. Ultimately, the 
construction is universal for such models and its essential feature, the lat- 
tice simplification, may point towards a renormalisation of the metric in 
continuum theories. 
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1 Introduction 



The application of tensor network methods has provided deep insight into topo- 
logicaUy ordered systems. Tensor networks provide a setting where the exotic 
characteristics of topological order (topology dependent ground level degener- 
acy, local indistinguishability of ground states, topological entanglement en- 
tropy, interplay with renormalisation, unusual realisation of symmetries) can be 
exactly studied. On the other hand, tensor network methods are sufficiently 
flexible for studying the excitations of these systems; in 2D, as is known, the ex- 
citations are quasiparticles with anyonic exchange statistics, which is the basis 
of Kitaev's topological quantum computer architecture [ ] . 

Among lattice models with topological order, Kitaev's quantum double mod- 
els [1] have a distinguished position. They were the first models in which the 
possibilities of the topological setting for quantum computation purposes were 
discussed, and include anyon models universal for quantum computation by 
braiding; on the other hand, quantum double models bear an intimate rela- 
tion to discrete gauge theories, exhibiting a rich group-theoretical and algebraic 
structure that pervades the study of their ground level and excitations, and 
in particular determines the properties of the anyons and their computational 
power. While they are not as general as string-net models [_'], the models that 
describe all doubled topological fixed points on the lattice, steps towards a re- 
formulation of the latter in the spirit of quantum double models have been taken 
in [■'-), 4, T)]. Recent progress in the perturbative production of quantum double 
codes via two-body interactions has been reported in [ ]. 

Tensor network methods were first applied to topological models on the lat- 
tice in [7], where an exact projected entangled-pair state (PEPS) representation 
of a ground state of the toric code was given. In recent years, a growing number 
of tensor network representations for topological states have been developed 
following the lead of [ ] (e.g. [8, 9, fO, II, 4]), which are notably exact for 
ground states of fixed-point Hamiltonians (quantum doubles, string-nets), and 
moreover lend themselves to deep analysis of their symmetries [1.2]. Tensor 
network methods for systems with anyonic quasiparticles, independent of any 
microscopic structure, have also been developed [ : ■!, I -']. 

The multi-scale entanglement renormalisation ansatz (MERA) is a partic- 
ular kind of tensor network structure, developed by G. Vidal [15, 16], which 
builds on the coarse-graining procedure typical for real-space renormalisation 
flows, but introduces layers of unitary operators between coarse-graining steps 
so as to reduce interblock entanglement, a procedure called entanglement renor- 
malisation (ER). In practice, ER led to tractable MERA representations and 
algorithms for ground states of critical systems in ID, traditionally hard for 
tensor network methods, and the numerical applications of the method now en- 
compass a wide class of lattice and condensed-matter models both in ID and 
in higher dimensions. ER was first applied to topological lattice models in [ ], 
where exact MERAs were given for ground states of quantum double models, 
both Abelian and non-Abelian. In [ ] ground states of string-net models were 
also written in MERA form. 

While understanding ground states of many-body quantum systems is cer- 
tainly of the utmost importance, it is a bonus for a method to be able to 
explore at least the low-lying excitation spectrum. The purpose of this paper is 
to show how the principles of the ER approach to quantum double models can 



3 



be extended to account for the whole spectrum in an exact way. The picture 
that emerges is somehow complementary to entanglement renormalisation: by 
giving up the coarse-graining, that is, using unitary tensors throughout, the 
ER method turns into a reorganisation of degrees of freedom which is essen- 
tially graph-thcorctical and proceeds by modifying the lattice by keeping the 
topological Hamiltonian the same at each step; this procedure yields finally a 
Hamiltonian consisting essentially of one-body terms, or an essentially classical 
one-dimensional spin chain model. 

Naturally, the structure of the tensor network with unitary operations defines 
a unitary quantum circuit, which in this case simplifies the structure of the 
Hamiltonian down to mostly one-body terms. In this sense, the construction can 
be regarded as an instance of the broad idea of quantum circuits diagonalising 
Hamiltonians, introduced in [17]. 

In this vein, we obtain an explicit (and geometrically appealing) construction 
of the well known mapping of the toric code to two classical Ising chains; but 
we also cover the non-Abelian cases. Remarkably, the construction is virtually 
universal for all quantum double models, and can be understood as a series 
of moves simplifying the structure of the lattice into a series of 'bubbles' and 
'spikes', each one of which talks only to one qudit; each quantum double model 
possesses a canonical set of operations transforming the model on one lattice into 
the same model on the next lattice. Since the lattice structure can be considered 
as the discretisation of a metric, this might open the door to speculations about 
a continuum counterpart (including, perhaps, dimensional reduction in going 
from the topological Hamiltonian in 2D to one-dimcnsional classical chains). We 
remark that the recent paper [lb] uses related ideas to find numerical methods 
applicable to lattice gauge theories. 

The structure of the paper is as follows. The Abelian case of the toric code 
is discussed in section 2, where we recall the basic elements of the construction 
of [<S], and proceed to a detailed presentation of the disentangling method. In 
section 3, the analogous construction is developed for general quantum double 
models; particular attention is given to the models on the torus, where the 
structure of the Hilbert space has to be taken into account carefully. Section 4 
contains a discussion and conclusions. In appendix A some algebraic background 
underlying quantum double models is provided, and appendix B concentrates 
on the simplest non-Abelian quantum double, that of the group S3. 
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2 Entanglement renormalisation of the toric code 



We begin by describing the construction in the simplest Abehan case, namely 
the toric code. This will allow us to make a connection with known results, 
namely the mapping to a pair of classical Ising chains. 

2.1 Elementary moves for the toric code 

We start from an arbitrary two-dimensional lattice A, with qubits sitting at the 
bonds. We do not specify the topology yet. The Hamiltonian, to begin with, 
has the familiar form of Kitaev's toric code [1]: 

H^ = -Y,Bp-Y,Ay , (1) 

P V 

with 

Bp = (g) , Ay = (g) X, (2) 

i£P jev 

mutually commuting stabilisers with support on plaquettes P and vertices V. 
Any code state \tp) (i.e., any ground state of Ha) satisfies the local conditions 
Bp\ijj) = It/i), AvIV") = for all P, V; whatever information it carries is purely 
topological, encoded in the eigenvalues of nonlocal operators. 

The entanglement renormalisation scheme proposed in [ ] uses two kinds of 
elementary operations introduced in [19], the plaquette move and the vertex 
move. In that work only code states were considered, and the result of applying 
elementary moves was to decouple qubits from the code in known states. More 
precisely, plaquette moves (P-moves) decouple a qubit q from the lattice by 
putting it in a product state |0) with the rest of the system: 

|code state; w)a ^-^ |code state; w)a' <8) |0)q , (3) 

where A' is the lattice resulting from the deletion of the edge containing q, 
therefore fusing the two plaquettes to which q belonged into one plaquette. Here 
w stands for any nonlocal information in the code state, which is preserved 
provided the plaquette move is purely a bulk operation, that is, it acts on a 
contractible region. Vertex moves (V-moves) effect a similar simplification, 

|code state; w)a |code state; w)a" ® \+)q , (4) 

with A" obtained by contracting the edge containing q into a vertex (subsuming 
the two vertices of the contracted edge). In both cases the topological infor- 
mation carried by the code state is unchanged, if the moves have support on a 
contractible region. 

The programme of entanglement renormalisation can be carried out explic- 
itly and exactly for the code states, as shown in [ ], using P-moves and V-moves 
organised into disentanglers and isometrics (in the latter case the decoupled 
qubits being discarded from the system). This gives rise to an exact MERA 
ansatz for the code states, whereby the top of the tensor network contains only 
topological degrees of freedom; for instance, in the case on the torus, these topo- 
logical degrees of freedom are two qubits on edges along noncontractible loops 
of T\ 
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Figure 1: Deformation of the lattice after a P-move. The P-move consists of 
CNOTs (red arrows). Qubit q "decouples" if in the ground state, but the complete 
picture is that it remains in the deformed lattice, with its own plaquette (a bubble 
containing the flux initially in the left subplaquette). To be consistent with the non- 
Abelian generalisation to be discussed later on, the bubble is attached to one of the 
vertices talking to q. 




Figure 2: Charges under a P-move. Redistribution of electric charges and 
magnetic fluxes rrii in the P-move of figure 1. Note that the redistribution would be 
different (swapping fluxes mo, mi in the deformed lattice) if the control qubits there 
belonged to the plaquette Pi, although for the ground states the effect is the same. 
Each deformed vertex contains the same electric charge as its undeformed predecessor. 
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Let us consider more in detail a P-move ?7piaq (figure 1, left). This consists 
of simultaneous CNOTs with target qubit q and control qubits in the remaining 
edges around one of the two plaquettes Pq, Pi containing q (say, Pq). Call P 
the plaquette obtained from Pq and Pi by deleting the bond for q. Only four 
terms in the Hamiltonian change: 

C/plaq(Sp, + Bp, + Av„ + Av,)Ul^^ = Z, + ZgBp + Ay^ + Ay^ , (5) 

where vertices Vb and Vi talk to qubit q, and yj, V-[ are the new vertices after 
deletion of the bond containing q. 

But this calls for a more general interpretation than just the coarse-graining 
(3) of the ground level. The first term in the rhs of (5) describes qubit q having 
eaten the flux in plaquette Pqj but not otherwise subject to vertex constraints; 
the second term contains all qubits in Pq and Pi, including q counted once, 
so the new Hamiltonian is not a toric code with q removed, but a toric code 
with a deformation of the lattice still containing q in a bubble^ i.e., in its own 
plaquette Pq contributing to the total flux (figure 1, right). Note that we draw 
Pq attached to one of the vertices of the bond previously containing q; this is 
needed for a consistent generalisation to non-Abelian cases. At any rate, q is 
not subject to any vertex constraint since it talks twice to this vertex V^', and 

X2 = 1. 

The effect of the P-move in terms of charges and fiuxes is drawn in figure 2. 
The electric charges remain at their vertices, while the bubble Pq of q contains 
the flux across Pq, and the complementary new plaquette P{ contains the flux 
previously in Pi. Note that while the effect on the ground level is symmetric 
with respect to the choice of control qubits, this is not the case when there are 
nontrivial fiuxes. 

Note also that the total charge and flux of the entire region affected by the 
local operation C/piaq are conserved, as they should be. 

Elementary V-moves are entirely dual to P-moves in the toric code (and, 
more generally, in any Abelian quantum double model). To be precise, the 
duality is established by going over to the dual lattice and performing a global 
Hadamard, thus swapping X and Z operators. 

Figure 3 shows the effect of such a move, J/vcrt- The operation consists of 
CNOTs with the control qubit q in the bond from vertex Vb to Vi (the one to 
'decouple' in state |-|-) if in a ground state) and targets those sharing with q 
one of the vertices, say Vq. The result is that the Hamiltonian is mapped to the 
Hamiltonian of a lattice where the bond of q has been contracted, identifying 
its two vertices into a new vertex V/, to which is attached a 'spike', i.e., a bond, 
containing g, whose other vertex is not connected to any other bonds. As 
shown in figure 4, the fluxes in the new plaquettes Pq and P[ are inherited from 
those in Pq and Pi, while the electric charge at vertex contains the charge 
previously at vertex Vq, and the vertex VI at the end of the spike contains the 
charge of Vi . 

2.2 Disentangling the Hamiltonian: a simple example 

With the set of elementary moves at hand, we can disentangle the Hamiltonian 
by simplifying the structure of the lattice. 

Think of a toric lattice, to be specific an x Ly square lattice Aq with 
periodic boundary conditions, on which a toric code TC(Aq) is defined. By 
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Figure 3: Deformation of the lattice after a V-move. As before, red arrows 
stand for CNOTs. Qubit q belongs, after the operation, to a spike (a bond with one 
of its vertices not connected to any other bond). 




using P-moves and V-moves as described in the previous section, the lattice gets 
deformed into a sequence of lattices A„, n — 1, 2, . . ., with the same number of 
vertices, edges, and plaquettes, and hence the same topology, but increasingly 
simpler connectivity, and keeping the lattice system a toric code TC(A„) at 
each step. States flow in this process just by redistributing the charges and 
fluxes. 

The lattice can be simplified to reach a particularly appealing form, to which 
corresponds a Hamiltonian consisting of three sectors: 

• A trivial sector corresponding to two qubits encoding the topological de- 
grees of freedom. 

• An Ising chain of L^Ly — 1 qubits encoding the structure of magnetic 
(plaquette) fluxes. 

• An Ising chain of L^Ly—l qubits encoding the structure of electric (vertex) 
charges. 

To understand this, note that the minimal toric square lattice consists of 
two links along two nonequivalent nontrivial loops on T^ , beginning and ending 
at a unique vertex, and defining a unique plaquette. In this minimal case, 
the Kitaev Hamiltonian is trivial and all computational basis states are Kitaev 
logical states, there being no room for excitations in the Hilbert space. 

On performing P-moves and V-moves, the topology of the lattice is left un- 
changed, but the number of nontrivial cycles along links in the lattice can be 
reduced to the minimal case, whereby two qubits carry all topological informa- 
tion. Information about excitations in the original lattice will appear in the rest 
of the code qubits. 

Let us illustrate the simplification in the case of a toric lattice with four 
plaquettes and four vertices. The steps are depicted in figure 5. We consider a 
basis state of the Hilbert space determined by magnetic fluxes to^ for plaquettes 
Pij and electric charges Cij at vertices v^j, taking values in {±1} (being the 
eigenvalues of the corresponding plaquette or vertex operator). Indices i, j 
run from 1 to 2 labelling the horizontal and vertical logical coordinates of the 
plaquettes and vertices. By overall charge neutrality, Yiij ™ij = +1 = Yiij ^ij- 
Topological information is encoded in two nonequivalent nontrivial loops, the 
eigenvalues of the corresponding logical Z operators being zjj, zy. After a series 
of P-moves and V-moves, the topological data are carried by the two qubits 
shown in dark in figure 5, which will be denoted as qv.H- Qubits carrying 
magnetic information will be referred to as qM,i, i — 1, 2, 3, and are those 
whose links enclose a solitary plaquette. Qubits carrying electric information 
are shown in light colour in figure 5, will be denoted as qE.i^ i — 1: 2, 3, and 
are those whose link has a free vertex. The final Hamiltonian is 

3 

Affinal = - ^ Z{qM, i) - Z{qM,l)Z{qM,2)Z{qM,3) 

1=1 

3 

- J2Xi(lE,^) - X{qE,l)X{qE,2)X{qE,3) ■ (6) 
i=l 

The precise redistribution of fluxes and charges will be discussed in general later, 
but for the moment let us remark that, as promised, there are no terms in -fffinai 
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Figure 5: Simplification of the lattice for a 2 x 2 toric code. The end lattice 
features two topologically nontrivial links and a series of bubbles and spikes, all at- 
tached to the same vertex. We have kept qubits in their places and only deformed the 
lattice links. Red arrows stand for CNOTs and point from control to target qubits. 
Dark qubits are those that hold a full qubit of topological information. Other colours 
are used for qubits that hold the flux (medium dark) or charge (light) of a plaquette 
or vertex of the original lattice. 
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associated with qv,H, and that the magnetic sector (terms affecting the qM,i, 
i.e., involving Zs) and the electric sector (terms affecting the qE,i), i-e., involv- 
ing Xs) are each intimately related to an Ising chain Hamiltonian. A cleaner 
illustration of the lattice structure associated with Hamiltonian (6) is given in 
figure 6 (the Abelian nature of the code allows for generous rearrangement of 
bubbles and spikes). 




Figure 6: A clearer picture of the disentangled 2x2 toric code. Light qubits 
carry electric information, medium dark qubits carry magnetic information, and dark 
qubits carry topological information. The Hamiltonian for this code is (6). 

The correspondence with the Ising chains can be made explicit by especially 
simple P-moves and V-moves acting on the configuration of figure 6 by direct 
application of stabiliser formalism results. These are shown in figure 7. Denote 
as C(< -i— c) the CNOT operation with control c and target t. Starting from (6), 
the corresponding operation is 

C{qE.2 qEA)C{qE.3 qE,2)C{qM.3 ^ qM,2)C{qM,2 ^ qMs) , (7) 

whose application yields the doubled Ising Hamiltonian: 

-fffinal = - ZiqM,i) - Z{qM,i)Z{qM,2) - Z{qM,2)Z[qM,3,) ~ Z{qM,3.) 

- X{qE, i) - X{qE^ l)X{qE, 2) - ^(te, 2)^(te, 3) " X{qE, 3) • (8) 

2.3 Disentangling the Hamiltonian: full structure 

As for the explicit moves and redistribution of charges and fluxes, necessary to 
understand the unravelling of the Hamiltonian in general, consider an L x T 
toric square lattice as shown in figure 8. We label vertices according to their 
coordinates as v^ j , where i runs from to L — 1 and j runs from to T — 1 . 
Plaquette p^ has vertex v^ as its lower left corner. Qubits in horizontal links 
are denoted qh^^-, with vertex v^ to their left, and qubits in vertical links as 
qVj , pointing upwards from vertex 
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Figure 7: From the disentangled toric code to Ising chains. The sequence of 
steps to bring the lattice explicitly to the form where the toric code Hamiltonian is the 
sum (8) of two noninteracting Ising chains. Red arrows, as usual, stand for CNOTs. 
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Figure 8: General toric code simplification: setting. The starting point for 
the discussion of the full structure of the procedure of simplification of the toric code. 
Here qubits are split into two groups according to whether they belong to a horizontal 
bond (qhj^) or a vertical bond (qv^^) in the original L x T toric lattice. Electric 
charges Bij = ±1 at the vertices and magnetic fluxes rriij = ±1 across the plaquettes 
are labelled as well by obvious row and column indices, i = 0, 1,...,!/ — 1 and j — 
0, 1, . . . , T — 1. Overall charge neutrality on the torus means Ylu "m-ij = 1 = Ylu ^'i- 
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Consider the basis of simultaneous eigenstates of the set of all plaquette 
operators (eigenvalues to^ with Y[ij rriij — 1) and vertex operators (eigenvalues 

Bij with riij = 1)' together with logical Z operators Zh = Ili^o^ ^qh, o 

(eigenvalue zh) and Zy — nJ=o^ ^qvo j (eigenvalue zy)- 

The process of simplification of the lattice will be described by the P-moves 
and V-moves applied to a basis state to bring the lattice to a form generalising 
figure 6, and the identification of the qubits carrying the eigenvalue of each 
vertex, plaquette, and logical original operator at that stage. 

There is a huge amount of freedom in the choice of the process. I will 
be describing one possible choice, not optimal in any sense, but allowing a 
systematic treatment. 

Step 1: Create bubbles along rows Figure 9. 

Operation Z^i = U'°^ with [/™* = Pi^^Ps.j . . . Pl-i,j. P- 

move Pi^j is defined as the application of CNOTs on target qubit qVj ^ with 
controls the rest of the qubits of plaquette pi_i^j. Explicitly, 

Pij = C(qv^j ^ qh,;_i_j+i)C(qv,;^^ ^ qVi_i_ j)C(qv,_j ^ qh,_i,j) . (9) 

Step 2: Create spikes along rows Figure 10. 

Operation ^2 = Oj^o' C/"" with U^" ^''^ - Vi_2,jVL_3,j ■ ■ • Vq.,. 
Each V-move here consists of a single CNOT, V^^ j — C(qhj_,_]^^^- •(— qh,^^), since 
the qubit in the affected bubble (qVj_|_i j) is acted upon twice, and hence un- 
changed. 

At this stage, the toric lattice has been simplified to a single column, to 
the vertices of which are attached a collection of bubbles and spikes. The next 
moves reproduce the strategy followed before for each row. 

Step 3: Create bubbles along the column Figure 11. 

Operation = Pq . . . P^f-_2J where each P-move consists now of L CNOTs, 

Pj = C(qhi_i^^+i ^ qvi_j-) . . . 

X C(qhi_i_j+i ^ qvL_i,j)C(qhi_i^^+i ^ qhi-i.j) ■ (10) 

Step 4: Create spikes along the column Figure 12. 

Operation Ui = Vy_2 • • • where each V-move consists now of L CNOTs, 

= C(qhi_2,j+i ^ qvo^j) . . .C(qho^^+i ^ qVo^j-)C(qVo^^-+i ^ qv^j) . (11) 

After this set of operations, the topological information is carried by qubits 
qhi_i (eigenvalue zh) and qvp (eigenvalue zy)- AH magnetic fluxes in the 
original plaquettes except for rriL^i^T-i sue now located in onc-qubit bubbles, 
according to the rule 

flux rriij is carried by qubit qVi^ij, < i < i — 2, Vj, 

flux rriL-i.j is carried by qubit qh^_]^ .^^j^, \fj. (12) 

Flux mL-i,T-i ensures overall flux neutrality and is carried by the large pla- 
quette on whose perimeter lie all qubits (with the caveat that qubits carrying 
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Figure 9: General toric code simplification, step 1. Creating bubbles along each 
row by P-moves from right to left. 
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Figure 10: General toric code simplification, step 2. Creating spikes along each 
row by V-moves from left to right. 
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Figure 11; General toric code simplification, step 3. Creating bubbles along 
the remaining column by P-moves from top to bottom. 
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Figure 12: General toric code simplification, step 4. Creating spikes along tlie 
column by V-moves from bottom to top. 
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topological information appear twice, so do not contribute to the associated 
stabiliser) . Similarly, all electric charges except for eo, o appear now in vertices 
at the ends of one-qubit spikes, according to the rule 

charge eoj is carried by qubit qvg ^ < j < T — 1, 

charge Cij is carried by qubit qh^_i ,j, 1 < i < i — 1, Vj, (13) 

and charge eo, o appears at a vertex shared by all qubits (twice by both qubits 
carrying topological information). The final lattice is drawn in figure 13. 




Figure 13: General toric code simplification, step 5. Final lattice. 

Hence, the operation diagonalises the Hamiltonian in the sense 

that (a) logical information does not appear in the Hamiltonian but is encoded in 
two physical qubits, (b) all independent eigenvalues of stabilisers in the Hamil- 
tonian appear now as eigenvalues of one-body operators. Global topological 
charge neutrality is ensured by two many-body terms. As explained before, a 
transformation into two copies of an Ising chain can subsequently be performed. 

As regards the efficiency of the procedure, the number of CNOTs to be ef- 
fected is 0{LT). If operations are parallelised as much as possible, the operation 
can be performed in 0{L + T) steps. 
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3 Quantum double models 



In this section we apply the ER-inspired scheme described in the Abehan case 
of the toric code to general quantum double models 0(0) based on non-Abelian 
groups [1]. While this construction holds exactly, with a more natural structure, 
in the class of quantum double models D(iJ) based on Hopf algebras [4], which 
contains and generalises the D{G) models, we settle for a discussion of the group 
case since the language is bound to be more familiar. 

3.1 Hamiltonian, charges 

The quantum double model based on a finite group G, or D(G) model, is defined 
on an 2D lattice A where quantum degrees of freedom of dimension \G\ sit on 
the edges. Orthonormal bases {\g)} are chosen for each oriented edge, labelled 
by group elements; the bases for the two orientations of any given edge are 
related by the inversion g i— g~^- The Hilbert space of an oriented edge may be 
identified with the group algebra CG of complex combinations of group elements. 
The Hamiltonian of the D(G) model is 

Ha^-J2Bp-J2Av , (14) 

P V 

where the first sum runs over plaquettes of A, the second over vertices of A, and 
Ay, Bp are mutually commuting projectors with support on the edges adjacent 
to the vertex and along the boundary of the plaquette, respectively. 

Plaquette projectors select configurations where the product of group ele- 
ments along the plaquette boundaries is the identity element of G; 

Bp \gi, 32, ■ • ■ , 5r) = 6e{gr ■ • -9251) \9i, 92, 9r) , (15) 

where gi, . ■ ■ , gr are group elements defining a computational basis state of the 
edges along the plaquette boundary, oriented and ordered along an anticlockwise 
circuit with arbitrary origin, as illustrated in figure 14. Vertex projectors have 
the form ^ 

Av \gi, g2, ■ ■ ■ , gs) = -r^ \kgi, kg2, kgs) , (16) 
1^1 fcGG 

where gi , . . . , gs are group elements defining a computational basis kct in the 
bonds adjacent to the vertex, oriented towards it, as in figure 14. 

Simultaneous +1 eigenstates of all Ay and Bp are ground states, and excited 
states break the ground state conditions Av\ip) = Bp\^Ij) = at some 
vertices or plaquettes. We can speak of particle-like excitations living at broken 
plaquettes and vertices, and it turns out that these particles are characterised by 
topological charges from an anyon model [20]. Topological charges are properties 
of regions that cannot be changed by physical operations within the region, and 
can be determined from measurements on its boundary (think of electric charge, 
measured by Gauss's law); more generally, a topological charge label can be 
assigned to any closed loop independently of whether it is the boundary of a 
region or not, which is crucial for the study of topologically nontrivial systems, 
as we will see in the next section. 

Topological charges in the quantum double models can be understood alge- 
braically: they are labelled by irreducible representations of Drinfel'd's quantum 
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Figure 14: Quantum double Hamiltonian. Conventions for the definition of 
plaquette and vertex projectors in equations (15) and (16). 



double T){G), a quasitriangular Hopf algebra which can be computed for each 
group G [21]. More details will be given later; for the moment let us remark 
that the distribution of topological charge on the lattice may be characterised 
locally by assigning charge labels to lattice sites, where, following Kitaev, a 
site is a pair of a plaquette and one of its neighbouring vertices. A charge has 
thus in general a magnetic part, associated with the plaquette, and an electric 
part, associated with the vertex; a charge with both nontrivial magnetic and 
electric parts is called a dyon. The part of Hamiltonian (14) acting on a site 
assigns energy —2 to the trivial charge (the unique charge with trivial electric 
and magnetic parts), —1 to purely magnetic or purely electric charges, and to 
dyons. 

3.2 The smallest torus 

To generalise the treatment of the toric code given in section 2, we first need to 
understand the topological degrees of freedom of the D(G) model on the torus, 
which are not covered by the previous local characterisation of topological charge 
distributions. 

In the toric code, the ground level has fourfold degeneracy; usually this 
is interpreted as the set of states of two logical qubits associated with the two 
noncontractible loops on T^ as in the left part of figure 15. The logical Pauli Z^or 
and .^vert measure magnetic quantum numbers rrihor, "^vcrt along the horizontal 
and vertical loops; this amounts to determining 'half the topological charge' for 
each loop. This defines a basis |TOhor, ™vGrt)hv of eigenstates in the ground level. 

Alternatively, one may choose as commuting set of observables the logical 
Pauli Z\ior along a horizontal loop of the direct lattice, and the logical Pauli 
Xhor along a horizontal loop of the dual lattice (centre of figure 15). In this 
case one has a basis |(7hoi)h labelled by the full topological charge label ghor = 
(rniior, Ghor) associatcd with the horizontal loop, including the eigenvalue mhor of 
the magnetic Zhor operator and the eigenvalue Chor of the electric Xhor operator. 

If we choose the corresponding logical Pauli operators along the vertical loop 
(figure 15, right), then we have another basis |(?vort)v labelled by a full topological 
charge. From the abstract properties of an anyon model [ ], we know that the 
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change of basis between |ghor)h and |(7vcrt)v is the unitary operator known as 
the topological iS-matrix: 

I Qvci-t ) V = ^ I 9hor ) h •S'gi^, , g^ert • (17) 
9hor 

This relation encodes a sort of uncertainty principle for topological labels on 
T^: if we determine the full label along the horizontal loop, the vertical label is 
'maximally spread', and vice versa. The states of the hv-basis work as sorts of 
coherent states where a compromise in the measurements of the horizontal and 
vertical labels is made so as to have mutually commuting observables. 



Zho 



ZhoT 



Zvcrt 
|"^hor, ™vcrt)hv 



9lior)li — \ {m.hoi, ehor))h 



^vcrt -^vcrt 
^vcrt)v — K^verti ^vert))v 



Figure 15: Logical operators in the toric code. Sets of commuting observables 
yielding complete topological information in the ground level. Left: Two logical Pauli 
Zs measuring 'half the topological charge' along each noncontractible loop on T^. 
Centre: A pair of logical Pauli Z and X measuring the whole topological charge along 
the horizontal loop. Right: The same for the vertical loop. In each case, the associated 
bases for the ground level are displayed. 

When considering non-Abelian models, we should not expect the hv-basis to 
survive. Indeed, in an Abelian D(G) model the number of topological sectors 
is |Gp, a perfect square, and one can successfully construct the |r7ihor, "^vcrt^hv 
states due to the separation of electric and magnetic degrees of freedom. In a 
non-Abelian D(G) model one has in general a number of charges different from 
a perfect square, for instance eight charges in the 0(83) model. 

However, the h-basis and the v-basis can be generalised for non-Abelian 
models. More than that, their existence stems from the modular nature of the 
underlying anyon model and is therefore a fundamental property. 

As in the Abelian quantum double model on the minimal toric lattice 

consisting of two edges contains the full basis of the ground level; in contrast 
with the Abelian case, the Hilbert space of these two qudits is not identical 
with the ground level, that is, there is room for excitations even in such a small 
system. Exactly for which excitations is ultimately a representation-theoretical 
problem. 

To wit, let |a, b)ch be the elements of the computational basis for the minimal 
torus, with a, b G (hgure 16). The action of the unique plaquette and vertex 
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Figure 16; Computational basis for the minimal torus in the D(G) model, 
projectors reads 



Clearly, if G is Abelian these are always the identity operator, so the Hilbert 
space of the two qudits is identical with the ground level. But for non-Abelian 
groups, the conditions Bp = jV'), Ay = IV") are nontrivial. We can 
specify a state by giving two sets of charge information: 

• A horizontal charge label ghor corresponding to one horizontal nontrivial 
loop, as in the previous discussion. 

• A bulk charge label gbuik corresponding to the charge sitting at the site 
defined by the unique plaquette and the unique vertex of the minimal 
square, together with internal degrees of freedom of the bulk charge as 
dictated by the microscopic model. 

Not all pairs (ghor, gbuik) are compatible: in particular, some charges never 
appear in the bulk, while some bulk charges are incompatible with some, but 
not all, of the loop charge labels. The condition for such a pair to be compatible 
is that the fusion of the two charges contains the horizontal charge: 



The reason for this is consistency of the charge measurement along the loop: if 
the loop is deformed by sweeping the whole bulk of the torus, the same charge 
label should be measured in the resulting loop, because of periodicity. Since 
the new loop is equivalent to the concatenation of the original loop and a loop 
enclosing the bulk charge, condition (19) follows (see figure 17). 

The discussion of these phenomena for T){G) models requires some algebraic 
machinery and we leave it for appendix A, where we attempt a pedagogical 
account; the case of the smallest non-Abelian group S3 is considered in appendix 
B. Here we just list the results: 

• Charges in a D(G) model are labelled by irreducible representations of 
Drinfel'd's quantum double algebra 0(0). These are given by pairs (C, a), 
where G are conjugacy classes of G, and a are irreducible representations 
of the centraliser group of C; more details appear in appendix A. 



Bp I a, b)ch 



S{ab, ba) \a, b)ch , 




(18) 



giior X gbulk ghor 



+ 



(19) 
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Figure 17: Consistency condition for loop and bulk charges. The loop in the 
right hand figure is, by periodicity, equivalent to the first horizontal loop, so the charge 
measured there has to be the same. Since the final loop is obtained by concatenating 
the original loop and the circuit measuring the bulk charge, the condition ghor x 9buik — 
Qhor + • • . follows. 



• The Hilbert space in the smallest torus consisting of two edges has an 
orthonormal basis: 

{\{C,a;k, z;,)buik;(C',a')hor); {C,a)<E,{C', a') ^ {C , a') ® . . .} , (20) 

where fc, Vi are degrees of freedom determining specific vectors in an irre- 
ducible representation space for (C, a). The projectors onto definite bulk 
and horizontal charge labels are given in equation (45) in appendix A. 

• From the point of view of the Hamiltonian of the D(G) model, configu- 
rations with vacuum bulk charge (Ce, 1+) (where the conjugacy class is 
Ce ~ {e}, and 1+ is the trivial representation of G) are ground states, 
with energy —2. The energy of states with pure magnetic charges (C, 1+) 
(where C 7^ Ce, and 1+ is the trivial representation of the centraliser of C) 
and that of pure electric charges (Ce, a) (where a 7^ 1+) in the bulk is — 1, 
due to the breakdown of the vertex and plaquette condition, respectively. 
Finally, the energy of states with a bulk charge (C, a) that breaks both 
the plaquette and vertex conditions (C 7^ Ce, a =/= 1-|_) is 0. 

3.3 Elementary moves 

The generalisation of the Abelian P-moves and V-moves to non-Abelian quan- 
tum double models was given in the Appendix of the preprint version of [8], in 
the case of groups; essentially, the CNOT operations of the Abelian case turn 
into controlled multiplication operators of the type \g, h) 1— >■ \g, gh) and their 
variants obtained by changes of orientation of the edges for control and target 
qudits, and these operators are composed to build the elementary moves. For 
models based on Hopf algebras, the elementary moves, also constructed from 
controlled multiplications, were defined in [-1] and give rise to a MERA for the 
ground states and the disentangling procedure of this paper in an analogous way; 
but as already mentioned, we concentrate on the case of groups here. We will 
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not give general expressions, but concentrate on simple examples from which 
the general procedure is straightforwardly inferred. 

A plaquette move or P-move has to be specified more carefully than in the 
Abelian case, since controlled multiplication operators acting on a target qudit 
no longer commute with each other. In general, we need specify both a plaque- 
tte and one of its edges as a target, or equivalently, one of the vertices along the 
boundary of the plaquette. A convention necessarily links the specific multipli- 
cation operators and the lattice transformation, including the orientation of the 
edges: ours is defined in figure 18, which we choose to associate with the site 
(Po, ^o). 




Figure 18: P-moves for a general D{G) model. Solid arrows stand for controlled 
left group multiplication operators in the form \g, h) M- \g, gh); these are to be ap- 
plied in the order defined by the dotted arrow, assuming the orientation of the edges 
involved as in the figure. All labels are group elements; therefore we show kets in the 
computational basis. We give the action of vertex representations of the group (gauge 
transformations) by showing how elements h, g are represented in the vertices Vo, Vi 
at the boundary of the target edge. After the P-move, the result is clearly a quantum 
double model in the new lattice that we draw, where the target edge talks twice to the 
lowermost vertex Vq and no longer has contact with the uppermost vertex V{. This 
we call a P-move associated with site (Po, Vo)- 

The vertex moves or V-moves are defined analogously, with the difference 
that one uses controlled multiplications with a single control and multiple tar- 
gets, that commute with each other, so one does not need to apply these oper- 
ations in a certain order. In the models based on Hopf algebras, this commuta- 
tivity property is lost and one has to take care of the order of operations as well 
[-1] (the property that makes the group case simpler is the cocommutativity of 
the Hopf structure in group algebras). Our conventions are defined in figure 19, 
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illustrating a V-move which we associate with the site (Vb, Pq). 




Figure 19: V-moves for a general D(G) model. Solid arrows, as before, stand 
for controlled multiplications \g, h) n- \g, gh); here they need not be applied in any 
particular order. The V-move in the computational basis, as well as gauge transfor- 
mations at the affected vertices, are shown in the bottom row. After the move, the 
control edge is the only one in contact with vertex Vq, and therefore the only edge 
feeling the gauge transformation by h. The plaquette in which the 'spike' is drawn is 
arbitrary from the point of view of the Hamiltonian; but to be definite we associate 
this V-move with site (Vo, Po)- 

3.4 Disentangling quantum double models 

For the purpose of analysing the quantum double Hamiltonian (14), the only 
relevant questions are those asked by the vertex and plaquette projectors Ay, 
Bp, whose eigenvalues we may consider a distribution of "binary charges" Cy, 
rup taking values in {0, 1}. As in the Abelian case, the original Hamiltonian 
gets transformed into the quantum double Hamiltonian in the deformed lattice, 
and the "binary charges" get reassigned as depicted in figures 2 and 4. The 
geometrical moves conserve the number of edges (degrees of freedom), vertices 
and plaquettes, and the operations on the state are unitaries preserving the 
energy. 

The disentangling procedure for general quantum double Hamiltonians (14) 
has thus the same graph-theoretical structure as the Abelian cases. The key is 
that the elementary moves simply redistribute the vertex and plaquette projec- 
tors along with the local deformations of the lattice, according to the general 
rule 

UpiV)HAU;^^)=HA'^,,, , (21) 
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where A is the original lattice, and Ap^^^ are the deformed lattices obtained 
after the P- and V-moves. 

If we work on a topologically nontrivial surface, there appears a difference 
with the Abelian models. Consider again the case of the torus: eventually the 
situation simplifies into a lattice where just two "topological" edges go around 
the two nontrivial cycles, and only one plaquette and one vertex talk (twice 
each) to them. In the Abelian case the topological edges do not contribute to 
the dynamics, essentially because the bulk charge in the quantum double model 
on the small torus consisting of just these two edges, as discussed previously, is 
always trivial. In a non- Abelian case, the small torus may contain bulk charges, 
as constrained by condition (19), and these contribute to the energy of the 
disentangled quantum double model. 

More precisely, let us consider the two simple cases analogous to the disen- 
tangling of the toric code into mostly sums of one-body operators (figure 6), 
and the mapping to the double Ising spin chain (figure 7). 




Figure 20: Disentangled quantum double. The lattice for the expression (22) 
of the Hamiltonian mostly as a sum of one-body terms. In a ket in the computa- 
tional basis, 'magnetic' qudits are in states \gi), 'electric qudits' are in states \kj), and 
topological qudits are in states \t-2)- 

In the first case, illustrated in figure 20, the qudits separate into three 
classes: r 'magnetic qudits', whose edges enclose a whole plaquette; s 'elec- 
tric' qudits, whose edges end at a vertex with no other edges attached to it; 
and two topological qudits. We write a ket in the computational basis as 
. . . , fci, . . . , fc^; ^1, according to the orientation conventions of fig- 
ure 20 (where r = s = 3), where the {gi}, {kj} and li, I2 correspond to the 
magnetic, electric, and topological qudits, respectively. The Hamiltonian reads 

r s 

H=-Y,B, - B^-Y^A,- A^, (22) 

4=1 3=1 
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where the action of the different operators in the computational basis is 

Bi \gi, . . . , gr; ki, . . .,ks; (i, £2) 

= ^e{gi) \gi, gr; ki, ...,ks; li, £2) , 
Bo \gi, ■ . ■ , gr; ki, . . .,ks; £1, £2) 

= 5,{g^^ . . . g;H^H^H2ti) Iffi, . . . , g.; fci, . . . , fc,; £2) , 
Aj \gi, gr; ki, .. .,ks; £1, £2) 

= j4| X! 1^1' • ■ • 7 ffr; fci, • ■ • , fcj?2"\ . . . , fc^; £1, £2) , 

' ' nGG 

^0 I51, gr; ki, .. .,ks; £1, £2) 

= 1^ ^ \ngin^^, . . . , ngrii"^; nki, . . . , nks; n£in~^, n^n^^) . (23) 

Each Ai, Bj is thus a one-body operator, and the terms Aq and Bq couple 
all qudits (except for the electric qubits in the case of Bq). This generalises 
Hamiltonian (G) for the toric code on the same lattice. Note that the precise 
location of the spikes is immaterial. 




Figure 21: Classical chains for quantum doubles. The lattice for the expression 
(24) of the Hamiltonian mostly as a sum of two-body terms. The notation is as in 
figure 20. 

In figure 21 we give our conventions for the lattice generalising the double- 
Ising-chain construction for the toric code. The Hamiltonian is now 

r— 1 s— 1 

H=-Bi - - Br^o - - - A,,o , (24) 

1=1 j=i 
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where the different terms are given by 



Bi \gi, gr; ki, .. .,ks; 



= 5e{gi) Iffi, . . . , ffr; fci, ■ • • , fcs; ^i, (-2) , 



Bi^i+i \gi, . . . , 5^; fci, . . . ,ks; h, £2) 



= 5e{gigi+i) \gi, . . . , 5r; fci, • • ■ , fc^; ^1, ^2) , 



Br,0 Iffl, ■ • ■ , 5r; fcl, • ■ • , fcs; ^1, (-2) 



= (5e(g-i^2-i£r'^2^i) Iffi, . . . , 3.; fci, . . . , fc,; ^i, ^2) , 



Ai |gi, . . . , gr; ki, . . .,ks; I2} 



" irT X!'^!' S"-' ^1*^ ^ ^2, ...,fcj, fcs; £1, £2) , 



\9i, ■ ■ ■ , 5r; fci, • ■ • ,fcs; ^1, -^2) 



= TT7| X!'^!' ^1' ■■■■.nkj, kj+in \ fc^; ^i, £2) , 



^s,o Iffi, gr; ki, .. .,ks; £1, £2) 

= ^|ngin"\ . . . , ngrn~^] ki, k^-i, nkg] n£in~^ , n£2n~^) . 



(25) 



Terms -Bi, are one-body, and terms Bi^i^i and Ajj^i are two-body (and 
define spin models which are essentially classical since all these terms commute) ; 
the projector i?^, couples one of the magnetic qudits with the topological edges 
(that this is not a many-body operator is one of the simplifications with respect 
to the general Hopf algebraic case); and A^^q couples the topological qudits 
with one of the electric edges and all of the magnetic edges. This generalises 
the double-Ising-chain mapping, with the terms Br, 0, ^s, representing the toric 
boundary conditions. The topological degeneracy is recovered easily: because 
of the one-body projectors Bi, Ai breaking the degeneracy of the two-body 
chains, the ground states of (24) feature product states in all the magnetic and 
electric qudits; and the remaining degrees of freedom are just the topological 
qudits, that is, the edges of the minimal torus studied in section 3.2, with the 
Hamiltonian —Bp — Ay defined by the operators in (18). Explicitly, the ground 
states of (24) are 



where £1, £2 are mutually commuting elements of G. 
3.5 The algebraic meaning 

The significance of the elementary moves goes beyond the mere fact that they 
intertwine the quantum double Hamiltonians in the two lattices. More funda- 
mentally, accompanying the geometrical evolution there is a topological charge 
redistribution more detailed than just the reshuffling of binary magnetic and 
electric labels. 




1 



\e, 
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Figure 22: The quantum double algebra. Conventions for the definition of pla- 
quette and vertex representations (27), (28). 



Consider again the problem of the local distribution of topological charges 
in an arbitrary state of a quantum double model. Instead of characterising this 
distribution of charges by a collection of charge labels at different sites (recall 
that a site is a pair of neighbouring plaquette and vertex, and that charge labels 
in general do not split cleanly into a magnetic and electric part), we can choose 
to endow the system with an algebra of local operators generalising projectors 
(15) and (16). Plaquette operators are representations of the algebra of functions 
/ from CG to the complex numbers, 

Bpyif) Iffi, 92, 9r) = /(ffr • ■■929l) \9i, 92, ■ ■ ■ , 9r) , (27) 

and vertex operators are representations of CG, defined by 

Av,p{k) \gi, 92, . . . , 9s) = Av{k) 52, • ■ • , 9s) = \kgi, kg2, kg^) , (28) 

in the conventions of figure 22. It can be shown that the Bpy{f) and Av^p{k) 
operators acting on a site define a reducible representation of D(G), decomposing 
into the different charge sectors for that site. Hence, we can simply say that the 
charge distribution is determined locally by the reaction of the lattice state to 
operators (27) and (28), for the different plaquettes and vertices. This has the 
advantage that we can work with plaquettes and vertices separately, without the 
complications due to the intricate nature of dyonic charges. Note that vertex 
operators (28) do not depend on the plaquette; this property is lost in the more 
general Hopf algebraic setting of [4]. 

The deepest property of the elementary moves is that they are intertwineis 
for representations Av^p{g) and Bpvif) along the flow of lattice deformations 
they define. (Ultimately, these moves realise morphisms in a category of repre- 
sentations.) The most complete approach to this deep property would entail a 
full study of the transformations of Kitaev's ribbon operators [ ], the building 
blocks of the representations Av{g) and Bp{f), under elementary moves, but 
for this paper we will just cover the transformation of the latter, and do so using 
figures. 

The general form of the intertwining property is 

Us, As (g) = As> (g) Us, , Us, Bs{f) ^ Bs>{f) Us, , (29) 
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where the move is either a plaquette or a vertex move associated with site Sq, 
and s' is the image of s after the lattice deformation. This obviously implies 
the intertwining property for the Hamiltonian, but is much more general. The 
transformation of the sites is subtle, though, and the notion of site as strict 
neighbouring pair has to be relaxed. Let us consider the examples that we 
already used to define the general moves. 

Any elementary move deforms one edge, which talks to two plaquettes Pq, Pi 
and two vertices Vq, Vi. In our convention, moves are associated with one site, 
say So = {Po, Vq). Consider the P-move studied in figure 18, and the transfor- 
mation of representations Ay, p{g) and Bpy{f) encoded in the geometric trans- 
formation of sites. In general, sites transform as s = (P, V) ^ s' = (P', V') 
according to the images of P and V under the deformation of the lattice (first 
row of figure 23). However, sites of the form s — [Po,V) with V ^ Vq have 
a more subtle transformation, as shown in the last three rows of figure 23. In 
those cases, s' — (Pq^ , V^'); where Pg^ ^° is the plaquette Pg conjugated by 
the path from V' to Vg' along the boundary of Pg. Operationally, this means 
that representations Ay get transformed into Ay, and Bp gets transformed 
into an operator whose argument is the conjugation of the group element along 
the boundary edge of Pg by the product of group elements along the path V'Vq 
(this conjugation is invisible for the operators entering the Hamiltonian). 

V-moves can be analysed in much the same way. Consider the move in figure 
19, associated with site sg = (Pq, Vg). Sites transform again as s = (P, V) i-> 
s' = (P', V) except for those having Vg as a vertex and P ^ Pq (note that sg 
itself transforms into (Po,V"g')). We have s = (P, Vq) ^ (P'^'o'^^i\ Vg'), that is, 
the plaquette is conjugated by a path from Vq to V{ (see figure 24). 
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Figure 23; Transformation of sites upon P-moves. The P-move of figure 18 is 
associated with site sq. Sites having the plaquette Pq in common with so have a subtle 
transformation involving a path along the boundary of this plaquette. 
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Figure 24: Transformation of sites under V-moves. The move in figure 19 
is associated with site sq. Sites having the vertex in common with so have again a 
nontrivial transformation. 
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4 Discussion 



In summary, we have extended the entanglement renormahsation approach for 
ground states of topological lattice systems to a unitary quantum circuit disen- 
tangling procedure that diagonalises the Hamiltonian. We give this construc- 
tion for Kitaev's quantum double models based on groups, on the understanding 
that the construction generalises directly to more general quantum double mod- 
els based on Hopf C*-algebras. The extension of the Abelian construction to 
the non- Abelian case is informed by the understanding of the structure of the 
Hilbert space and energy levels in the toric topology, and we give an account of 
this with details and a specific example to be found in the appendices. 

While this ER picture extended to the full spectrum of the quantum doubles 
is geometrically pleasant, it is probably not a practical way to construct the 
relevant codes from a small instance by adding degrees of freedom, since it 
requires very nonlocal operations. 

Note that the double Ising chain lurking behind the toric code has been ex- 
posed by different authors. In particular, a duality mapping was given in [22] 
showing that the partition function of the toric code is identical to that of two 
noninteracting Ising chains. The construction presented here is graphic and ex- 
plicit, and lends itself to generalisation for non- Abelian quantum double models. 
It also uncovers an intriguing geometrical picture of renormahsation, proceed- 
ing through deformations simplifying the original lattice, that can probably be 
extended to general string- net models building upon the work in [' ]. 

On the other hand, what we achieve in the end is (up to boundary condi- 
tions) the diagonalisation of a Hamiltonian consisting of mutually commuting 
projectors. The approach of Ref. [ ] is a general method for diagonalising 
strongly correlated Hamiltonians in terms of a quantum circuit, and the present 
work can be regarded as an illustration thereof (the tensor network behind this 
construction being a unitary quantum circuit obtained from unitarisation of the 
ER scheme of [ ]). 

As remarked in the introduction, a lattice is essentially a discretisation of 
the metric in a continuum theory. One of the notable points of this construction 
is that the flow of the quantum circuit is universal for all the quantum double 
models, being encoded in the simplifying moves acting on the lattice structure: 
whether this points to a meaningful procedure in the continuum limit is a sug- 
gestive question. 

Work using related geometrical ideas, leading to numerical methods for the 
Z2 gauge theory, was reported in [ ] . Similar ideas are being used by R. Konig, 
F. Verstraete, and collaborators, working towards a more general renormahsa- 
tion scheme for lattice systems, to be found in [-'■>]. 
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A Topological charges in D(G) models 



In a D(G) model, topological charges are given in terms of the group structure 
and its representation theory (in the recent paper [24], it was shown that not all 
such charges are physically distinguishable; however, we will not be concerned 
with these issues here). Charge labels are pairs (C, a), where: 

• C runs over conjugacy classes of the group, that is, orbits of group elements 
under the adjoint action t gtg^^. 

• a runs over irreducible representations (irreps) of the centraliser group 
Nc of the conjugacy class C. 

The centraliser Nc is constructed as follows: given c G C, its centraliser Nc 
is the subgroup of G consisting of all the elements of G commuting with c, 

Nc = {ne G: ncn,-^ = c} . (30) 

It is immediately checked that 

neNc ^ gng-^ G Ng^g-i , (31) 

so the centralisers of elements of the same conjugacy class are conjugate (and 
hence, isomorphic) to each other, justifying the definition of the abstract Nq- 
The size of this group is |iVc| = |G|/|C|. 

Pairs (C, a) correspond to irreducible representations of the quantum double 
D(G); for the purposes of this work, it suffices to present D(G) as the tensor 
product (g) CG of the space of complex functions on the group and the space 
of complex linear combinations of group elements, or group algebra. This vector 
space is spanned by elements Phg = 5h® where g, h G G, and 6h '■ k ^ 5h,k 
are Kronecker delta functions. The algebra structure in D(G) is given by the 
multiplication 

Phg ■ Pki ^ Shigkg-^) Phge . (32) 

An irreducible representation (C, a) of D(G) realises elements Phg as endo- 
morphisms of a linear representation space this can be identified with 

CG V", where V" is the representation space for group representation a. A 
basis of yC^'") is given by {|fc; vt)}, with k G C and {\vi)} a basis of V"', such 
that 

D'^{n)\v.,)^Y.\'^,)Df,in) , n e Nc . (33) 
j 

The matrix form of (G, a) reads [25, 26] 

D^'''''HPk9) \k, V.) = Shigkg-') Y.\gkg-\ v,) Df.igl^J , (34) 

j 

or, in other words, 

D^MiPhg) = Mgkg-') Shigkg-') D^MnJ ■ (35) 

Here fco is an arbitrarily chosen reference element in conjugacy class G, and a 
is chosen as a representation of the normaliser Nkg. The notation glN^g means 
a projection of g G G onto N^g, defined as follows: given the reference fco € G, 
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define a system of group elements x\^^ G G, one for each k £ C, such that 
a;^^fco(a;^^)~^ = k, and in particular x'^° = e. It can be checked that x'^^Nkg are 
mutually disjoint cosets whose union is G. Now the adjoint action of g maps k 
to gkg~^, and from this follows that glN^g = i^kt^ )~^9^ka element of 

Nkg, which is what appears in (34) and (35). Different choices of ko and 
give rise to equivalent representations. 

The left regular representation, that is, the representation of D(G) on D(G) 
defined by left multiplication 

a i-> L{a) , L{a) -.h^ab , (36) 

is reducible and contains each one of the irreps (G, a) of D(G) with a multiplicity 
equal to its dimension |G||a|, just as for finite groups: 

D(G)« K^^'"). (37) 

(C, q) a=l 

Explicitly, the following elements of D(G): 



[(C, a) fc, ^; = MM D'^,{{xlX'n-'xl)nxl{xiX\ (38) 

where k, m G C, are a basis of the quantum double algebra satisfying 

Phg ■ [(G, a) fc, £; z, s] 

= 5h{gkg-') ^[(G, a)gkg-\ i- j, s\ D%{{xff'Y^gx\g) , (39) 

and hence, they realise decomposition (37) if we identify index a with the pair 
(£, s). 

Projection onto the different (G, a) sectors in D(G) is achieved by multipli- 
cation with the following central algebra elements: 

2(C,.)^Mj^ ^ E x"(«)-^--^<)P.n, (40) 

where x" is the character of representation a. It can be checked that they form 
a decomposition of the unit into a sum of projectors: 

g(C, a) . g(C', a') ^ ^(C, a), (C, «') g(C, a) ^ ^ g(C, a) ^ ^^^^^ = J] P/.e . 

(C, ct) /i 

(41) 

Upon multiplication, they project onto the irrep blocks 0^ Va'^' of decompo- 
sition (37): 

Q(^' • [(G', a') /fc, £; i, s] = 6^^' ^^''"'^ [{C , a') k, z, s] . (42) 

The key algebraic structure defined on the D(G) lattice models is the algebra 
of ribbon operators [ ] (see also [ ] for a very readable introduction). Without 
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going into the details of the definition, these operators are assembled from one- 
edge operators L±{g) representing G (or, more generally, the group algebra 
CG), and T±{f) representing the dual algebra C*^ of complex functions on G: 

L+{g)a^ga, L^{g)a = ag^^ , T+{f) a ^ f{a) a , T^{f) a ^ f{a^'^) a . 

(43) 

The L± are just the left and right regular representations of CG, while the T± 

are equivalent to the left and right regular representations of C*^ (through the 

map g \/\G\ 5g-i , which is an electric-magnetic duality [•"]). 

By forming combinations of these elementary operations using the properties 

of CG and C*^ as Hopf algebras, one can build representations of D(G) and 

its dual. In particular, one can define projectors to 'measure the topological 

charge' associated with geometrical closed ribbons by representing the elements 
g(c,a) g 

Coming back to our smallest torus, we define topological labels for a horizon- 
tal loop, a vertical loop, and the bulk site, using the following representations 
of D(G): 

FhoriPhg) |a, b)ch = 5h{gag^^) \gag^^ , gb)ch , 
F^citiPhg) \a, b)ch = Sh-i{gbg^^) \ga, gbg''^)ch , 
FbuikiPhg) |a, fo)cb = 6higb-^a-^bag-^) \gag'\ .96g"^)cb • (44) 
From these representations, the projectors onto definite topological sectors are 

Fhor(S(^'"') \a, bU = W E ^''W E x'^Hj-'n-'xl) \a, nb),^ , 



|G| 
|G| 



keC neNk 



Fv„.t(Q(^'"0 \a, bU - W E ^'c-i(6) E x"{{4o)-'n-'xl) \na, &), 



keC n£Nk 



i^buik(Q('^'"')|a, b) 



J2S,{b-'a-'ba) x"(«)-^n-ix^J|nan-i,n6n-i) 



IGI 

' ' fcec neNk 



cb 



(45) 

Projectors Fbuik(Q'-*^' commute both with horizontal and with vertical 
projectors in (45). This justifies a labelling of states in the small torus Hilbert 
space by a bulk charge and a loop charge (horizontal or vertical). The bulk 
charge is relevant for the energy, since the Hamiltonian penalises the departure 
from trivial charge in the bulk; loop labels, on the other hand, do not play a 
role for the Hamiltonian. 

The ground level of the quantum double Hamiltonian is the -1-1 eigenspace 
of the projector 

^^buik(Q<^^^ \a, &)cb = 77^ E Se{b-^a-^ba) \nan-\ nbn-^),^ . (46) 

' ' neG 

It is easy to check that the dimension of this vacuum subspace equals the number 
of topological charge labels. We can give an orthonormal basis labelled by a 
charge label in the horizontal loop, i. e., the h-basis, as 

|(G, a))h = ^ E E xH«r'^''0 \k, ri)cb , (47) 
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or the alternative v-basis labelled by vertical charges, 



a)) 



1 




) \n, k ^)cb ■ 



(48) 



n 



The overlap of basis states of (47) and (48) gives the topological S'-matrix 
of the theory: 



where fep, io are the reference elements of conjugacy classes C, C", respectively. 

For the whole Hilbert space of the small torus one may expect to find a 
basis labelled by a loop charge, say the horizontal label, and a bulk charge (plus 
possibly internal degrees of freedom for the latter) . This is a good strategy, but 
in general not all such charge pairs are compatible. For Abelian G only the 
trivial charge may appear in the bulk. For non- Abelian models, some charges 
can never appear in the bulk, and additionally some bulk charges are excluded 
depending on the loop charge label. Exactly which charges may appear in the 
bulk is a property of the anyon model and ultimately a representation-theoretical 
problem. As argued in the text, the compatible pairs of bulk charge and loop 
charge are those that can fuse into the loop charge. If the group G is Abelian, 
the representation theory of T)(G) is Abelian and charges are only invariant if 
they fuse with the vacuum charge (Ce, 1+); this singles out the vacuum label 
as the only possible bulk charge^. For non- Abelian models, the allowed pairs 
must be computed from the fusion rules. We give the full details for G = S3 in 
appendix B. 

A comment on the notion of charge in the bulk is in order. The projectors 
^buik(S^*"'"'') probe the topological charge in a contractible region, namely the 
square defined by the edges (square aba~^b~^ in figure 25); this happens to cover 
the whole area of the torus but cannot see topological quantities, because it does 
not include the boundary conditions. If one measures nontrivial topological 
charge on the bulk, it is exactly compensated by an opposite charge sitting at 
the boundary, and measured along the opposite circuit {bab~^a~^ in the notation 
of figure 25). All in all, the total charge on the whole torus corresponds to a 
measurement along a circuit equivalent to that in the right part of figure 25, that 
is, a loop always homotopic to a point, which always yields the trivial charge. 

The way in which bulk charge appears in the non- Abelian D(G) models is 
one aspect of Kitaev's construction which differs from a theory of gauge fields. 
In a discrete gauge theory, edge variables are gauge variant; therefore, toric 
boundary conditions only fix the periodicity of the edge variables up to a gauge 
transformation defined on the boundary of the bulk, which in turn is related to 
the bulk charge (in, for instance, a U{1) gauge theory, single- valuedness of the 
gauge transformation imposes Dirac's quantisation condition, which constrains 
the value of the bulk magnetic fiux). In Kitaev's lattice models, however, edge 

^More consequences can be drawn; for instance, since the magnetic part of the bulk charge 
has to be a conjugacy class of the commutator subgroup G' = {b~^a~^ba; a, b £ G), any 
charge not satisfying that condition never leaves charges invariant upon fusion; but it is not 
the point of this paper to develop the theory of fusion further. 



S(C',a'),{C,a) =v((C", a')\{C, a))h 



= u^E E x"'(«)-^^4)x"(«)-^^<) 



(49) 



kec eecnNk 
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Figure 25: Definition of bullc charges on the smallest torus. Cutting a torus 
into a square defines a contractible region, with a boundary along the circuit aba'^b^^ . 
The charge measured in this region is what we call bulk charge. The question of the 
total charge on the surface of the torus corresponds to the circuit aa~^bb~^, yielding 
always the trivial charge. 

variables correspond to physical states of physical qudits, which are strictly 
periodic on the torus. 
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B S3 and the 0(83) model 



In this appendix we illustrate the general formalism for the smallest non-Abelian 
group and the associated quantum double model. 

Let G = S3 be the group of permutations of three objects. It has six elements, 
organised into three conjugacy classes: the unit e; the transpositions to — (12), 
h = (23), <2 = (31); and the 3-cycles c+ = (123), c_ = (132). 

The topological sectors of the 0(83) model, that is, the irreducible represen- 
tations of the quantum double algebra 0(83), are labelled by pairs (C, a), where 
C is one of the conjugacy classes Ce = {e}, Ct = {to, ti, ^2}, and Cc = {c+, c_}, 
and a is an irrep of the centraliser Nc of the conjugacy class C. Since the cen- 
tralisers of the three conjugacy classes of S3 (shown in table 1) are Nc^ = S3 
(with irreps 1_|-, 1_, 2), = (with irreps 1+, 1_), and Nc^ = Z3 (with 
irreps 1+, 1(^, l^j), we have eight charges. 



Nc. - S3 



to 



to 



c+ 



1 



1 



1_ 



1 

1 



-1 



-1 



-1 



1 



1 

1 



w 

Lb 



u) 
w 



w 
u) 



Nc, = Z2 


e 


a 


1+ 


1 


1 


1_ 


1 


-1 



Nc. = Z3 


e 


b 


b^ 


1+ 


1 


1 


1 




1 




U) 


1^ 


1 


Q 


UJ 



Table 1: Irreducible representations of the centralisers 83, Z2, and Z3 of the 
conjugacy classes of 83. Here a; = 



The fusion rules of these representations (i. e., the Clebsch-Gordan decompo- 
sitions of their tensor products) are best presented in the manner of Overbosch 
[2()]. Labelling irreps by their dimension, we write 1 = (Ce, 1+); 1 = (Ce, 1-); 
2^,v,^,n' = (Ce, 2), (Cc, l+,u>,u>); and 3"''' = (C*, 1±). Then the nontrivial fusion 
rules can be written as 

1®1^1, 1(g) 2"= ^2^, l«)3"->3^ 
2^(8)2"^ ^ lei® 2^, 2^(g22'^2^®2'", 2^ (g) 3° ^ 3" ® 3^ 
3°(g)3'' ^ 1®2^®2''®2^®2"', 3" «) 3^ ^ 1 ® 2^ ® 22' ® 2^ ® 2^", (50) 

and the rules obtained by permutations in a, b and in x, y, z, w. In particular, 
each charge is its own conjugate. 

In the following we use the shorthand (e±) = (Ce, 1±); (e2) = (Ce, 2); 
{t±.) — {Ct, 1±), and {erf) = {Cc, 1?;) for the irreps of 0(83), where r] is one of 
the cubic roots of unity 1, uj, ui. 

For the explicit construction of these irreps we choose reference elements 
e & Ce, to & Ct, and c+ € Cc- The nontrivial fixed elements x\^^ such that 
x^^fco(x^^)~^ = k are chosen as = c+, a;^^ = c„, Xc^ = to- Projectors (40) 
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onto irrep blocks of the left regular representation read 



Q(e±) 


= i Pe (e ± to ± tl ± *2 



+ C+ +c_) , 




Q(e2) 


= iPe(2e-c+-c_) , 






Qit±) 


= i {Pt„ (e ± to) + Pt, ( 


e±ti) + Pt, (e±t2)} , 






= ^ {-Pc+ (e + r?c+ + ?7c. 


-) + Pc. (e + r;c+ + ??c_)} . 


(51) 



In the quantum double model defined on the smallest torus of figure 16, there 
are eight simultaneous +1 eigenstates of Ay and Bp (or just +1 eigenstates of 
the projector onto vacuum bulk charge, AyBp — Pbuik(Q''^'^'')): 

|e, e)cb, {|e, to)cb + |e, <i)cb + |e, t2)cb}, { |e, c+)cb + |e, c_)cb}, 

{\to, e)cb + \tl, e)cb + 1*2, e)cb}, {|*0, *o)cb + 1*1, *l)cb + |*2, *2)cb}, 

{|c+, e)cb + |c-, e)cb}, {|c+, c+)cb + |c-, c_)cb}, {|c+, c_)cb + |c-, c+)cb} ■ 

(52) 

As expected, there are as many as there are topological charge labels. 

We can rearrange the vacuum states (52) into the h-basis labelled by topo- 
logical charges associated with the horizontal loop in the following way: 

|(e±))h = {|e, e)cb ± |e, *o)cb ± |e, *i)cb ± |e, *2)cb + |e, c+)cb + |e, c_)cb}, 

I(e2))h = {2|e, e)cb - |e, c+)cb - |e, c_)cb}, 

|(<±))h = {1*0, e)cb ± 1*0, *o)cb + 1*1, e)cb ± 1*1, *l)cb + 1*2, e)cb ± |*2, *2)cb}, 

l(c?7))h = {|c+, e)cb + ^|c+, c+)cb + '7|c+, c_)cb + |c-, e)cb + ?7|c-, c+)cb + 7?|c-, c_)cb} • 

(53) 

These are exactly the states (47). The construction of the v-basis (48) is anal- 
ogous, and the overlaps of the elements of both bases determines the ^-matrix 
of the 0(83) model [2(i], which we list in table 2. 

The 36-dimensional two-qudit Hilbert space contains thus eight ground states 
and 28 excited states. The latter can be characterised as breaking the plaquette 
condition, the vertex condition, or both, yielding a magnetic charge, an electric 
charge, or a dyon at the single site defined by P and V. These states we char- 
acterise as the -1-1 eigenvectors of the projectors Pbuik(Q^'"' "^)- Let us study 
the excited states according to their bulk charge: 

• Bulk charge (e— ), that is, pure electric charges of type 1_: four states 
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X, 






6 X S(^c',a'),{C,a) 
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li) 


ID 
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1 


1 — 1 
u 






v((e+)| 


1 


1 


2 


3 


3 


2 


2 


2 


v((e-)| 


1 


1 


2 


-3 


-3 


2 


2 


2 


v((^-^2)| 


2 


2 


4 








-2 


-2 


-2 


v((/+)l 


3 


-3 





3 


-3 











v((t-)l 


3 


-3 





-3 


3 











v((cl)| 


2 


2 


-2 








4 


-2 


-2 


v(M| 


2 


2 


-2 








-2 


-2 


4 


v(M| 


2 


2 


-2 








-2 


4 


-2 



Table 2: Topological S'-matrix for the D(S3) model, computed from the overlap 
of the h-basis and the v-basis. 



with horizontal charge labels (e2), and all of the {cr}): 

|(e-)buik; (e2)hor) = {|e, c+)cb - |e, c_)cb}, 
|(e-)buik; (c?7)hor) 

= {|c+, e)cb - |c-, e)cb + j?|c+, c+)cb - ??|c-, c_)cb + ?7|c+, c_)cb - ??|c-, c+) 

(54) 

Notice how the horizontal charges are those which can remain unchanged 
upon fusion with the bulk charge. 

• Charge (e2) (pure electric charge of type 2): six states organised into three 
two-dimensional irreducible spaces: 

|(e2 : ?;i)buik; (e2)hor> = {|e, to)cb + w|e, ti)ch + w|e, i2)cb}, 

|(e2 : U2)buik; (e2)hor) = {|e, io)cb + w|e, ti)cb + w|e, t2)cb}, 

|(e2 : ui)buik; (i±)hor) 

= { |io, e)cb + e)cb + w|f2, e)cb ± |io, io)cb ± ti)cb ± w|i2, i2)cb}, 

|(e2 : i;2)buik; (i±)hor) 

= {1^0, e)cb + e)cb + w|f2, e)cb ± I^O, to)cb ± tx)ch ± ^1*2, i2)cb} • 

(55) 

Again the only horizontal labels arc the irreps that can fuse with (e2) 
without change, namely (e2) and (i±). 

• The charges (t±) never appear in the bulk. The transpositions to, t\, t2 
do not belong to the commutator group of S3; on the other hand, these 
representations always change the charges that they fuse with. 
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• Each charge (cl), (ccj), (clj), whose magnetic part is a flux of 3-cycle type, 
appears in the bulk in six states, organised into three two-dimensional 
modules. Let ij, as usual, run over the cubic roots of unity 1, w, co. Then 



{1^0, ii)cb + vlh, t2)ch + vlh, to)ch ± 1^0, c+)cb ± vlti^ c+)cb ± vlh, c+)cb}, 



I (cry : c+)buik; (i±)hor) 
1 

|(c77 : c_)buik; {t±) nor/ 

= {\to, t2)ch + vlh, to)ch + fj\t2, il)cb ± \to, C-)cb ± vlh, C_)cb ± V\t2, C_)cb}, 

|(C?7 : C+)bulk; (c77)hor) = {|c-, io)cb + ^|c-, tl)cb + t2)ch}, 



|(c?7 : c_)buik; (cf7)hor) = {|c+, io)cb + ??|c+, ii)cb + ??|c+, i2)cb} 



(56) 



Indeed, the only charges that can fuse with (cq) without change are (cry) 
itself and the (i±). 

This label distribution is summarised in table 3. 
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Table 3: Bulk and horizontal loop charge label pairs realised on the small torus. 
The entries represent the number of states belonging to the corresponding charge 
pair. 



For each bulk charge, an analogue of the topological S'-matrix can be con- 
structed by computing the overlaps of elements of the horizontal and vertical 
bases: the results are gathered in table 4. 
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Table 4: Analogues of the topological 5-matrix for the bulk excited states in 
the 0(83) model. By adjusting the phases in the definition of the bases in the 
last table we can obtain exactly the coefficients of the second table (all two- 
dimensional irreps of D(S3) have the same behaviour with respect to modular 
transformations) . 
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